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ABSTRACT 

The recent detection by the Fermi 7-ray space telescope of high-energy 7-rays from 
the radio galaxy NGC 1275 makes the observation of the very high energy (VHE: 
E > 100 GeV) part of its broadband spectrum particularly interesting, especially for 
the understanding of active galactic nuclei (AGN) with misaligned multi-structured 
jets. The radio galaxy NGC 1275 was recently observed by VERITAS at energies above 
100 GeV for about 8 hours. No VHE 7-ray emission was detected by VERITAS from 
NGC 1275. A 99% confidence level upper limit of 2.1% of the Crab Nebula flux level is 
obtained at the decorrelation energy of approximately 340 GeV, corresponding to 19% 
of the power-law extrapolation of the Fermi Large Area Telescope (LAT) result. 

Subject headings: galaxies: Seyfert — galaxies: individual (NGC 1275, 3C 84, Perseus 
A) — gamma rays: observations 



1. Introduction 



The search for 7-rays from radio galaxies is important for the understanding of the dynam- 
ics and structure of jets in active galactic nuclei (AGN). Radio galaxies are the only non-blazar 
extragalactic objects detected in the very high energy range so far. Even though radio galaxies 
contain jets and are considered AGN, their jet is not oriented toward the observer and therefore 
the radiation produced by the jet is not Doppler-boosted towards higher energies and luminosities, 
making them more challenging to detect in the very high energy (VHE: E > 100 Ge V) regime. The 



disco very of VHE 7-rays from the radio galaxy M 87 b y the HEGRA collaboration (jGotting et al 



200J), detected later by VERITAS (lAcciari et alJl2008af >. and from NGC 5128 (Centaurus A) by the 
HESS collaboration (j Aharonianl 12009) has shown that non-blazar AGN can produce very energetic 
photons from non-thermal processes. 

Radio galaxies are classi fied into two main families based on the morphology of their radio emis- 
sion (jFanaroff &i Rileylll974l ). whether it is core dominated (FR I) or lobe dominated (FR II). The 
differences bet ween these two familie s can be found in the radio energetics and in the discrete spec- 
tral properties (jZirbel &: Baumlll995l ). The large number of features that FR I radio galaxies share 
with BL Lac type blazars suggests a possible unification between the two sub-c lasses of AGN, in 
which FR I radio galaxies are BL Lac objects observed at larger jet viewing angles (jUrrv fc Padovani 
199Si ). 



Evidence for synchrotron emission in radio to X-ray energies from both the extende d structures 
and the core is well explained by relativistic particles moving in a beamed relativistic jet (jGhisellini et al 
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19931 ). A commonly considered mechanism for HE-VHE ( HE: high energy, 1 00 MeV< E <100 GeV) 
radiation is the synchrotron-self-Compton (SSC) process (| Jones et al.lll974l ). where the optical and 
UV synchrotron photons are up-scattered by the same relativistic electrons in the jet. Predictions 
conc erning the inverse Comp ton (IC) component have long been established for the 7- ray emis- 
sion (jBloom fe Marscherlll996l ) and frequency-dependent variability (jGhisellini et al. I ll989h . Besides 
leptonic scenarios, several models also consider a hadronic origin for non-thermal emis sion in jets. 
Acce lerated protons can initiate electromagne tic cascades or photomeson processe s (IMannheim 



1993), or directly emit synchrotron radiati on (lA haronian 



2002; 



Reimer et al 



2004) and pro duce 



7-rays through collisions with ambient gas (|Beall & Bednareklll999l : IPohl & Schlickeiserl I2OO0} ) . 



Modelling the blazar jet emission with a homogeneous SSC mechanism may imply particularly 
high Lorentz factors, T > 50, with consequent high Doppler factors and small beaming angles 9 ~ 1° 
(jKrawczynski et al.l 12002). Such a small beaming angle is in conflict with the unification scheme 
according to which FR I radio galaxies and BL Lac objects are the same kind of object observed 
at different viewing angles. Moreover, these high values for the Doppler factor are in disagreement 
with the small appare nt velocities ob served in the sub-parsec region of the TeV BL Lac objects 
Mrk 421 and Mrk 501 (|Marscherlll999l ). These considerations suggest a more complicated geometry, 
for example a decelerating flow in the jet wi th a consequent gradient in the L orentz factor of the 
accelerated particles and a smaller average V (jGeorganopoulos fe Kazanas!l2003l ) . As a result of this 
gradient, the fast upstream particles interact with the downstream seed photons with an amplified 
energy density, because of the Doppler boost due to the relative Lorentz factor r rc i. The IC process 
then requires less extreme values for the Lorentz factor and allows larger values for the beaming 
angle. In a similar way, a jet spine-sheath structure consisting of a faster internal spine surrounded 
by a slower layer has been also s uggested for the broadband non-thermal emission of VHE BL 
Lac objects (jGhisellini et al. 20 05]) . This mo del is supported by radio maps which show a limb- 
brightened morphology (jGiroletti et al.ll2004l ) and may explain the HE-VHE emission observed in 
radio galaxies at larger angles (#i aye r = I/Flayer ~ 20°). Observation of the VHE component from 
radio galaxies is therefore significant for the multi-structured jet modeling, as it can be related to 
the external lower structure of the jet itself. 

NGC 1275 (Perseus A, 3C 84) is a radio galaxy located in the center of the Perseus cluster 
and is one of the most unusual early-type galaxies in the nearby universe (z = 0.018). Its radio 
emission is core dominated, but it also has strong emission lines. In addition, the emission line 
system shows a double structure, corresponding to both a high-velocity and a low-velocity system. 
The puzzling nature of NGC 1275 makes it difficult to definitively classify it in a standard AGN sub- 
class. Its spectral energy distribution (SED) extends from radio to HE 7-rays. A source coincident 
with NGC 1275 at high con fidence was recently detected in HE 7-rays by the Fermi 7-ray space 
telescope (jAbdo et al.ll2009al ) . The time-average spectrum observed between August and December 
2008 is described by a power law in the energy range from 100 MeV to 25 GeV : 
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— = (2.45 ± 0.26) x 1(T 9 
dE 



E \ 



-2.17±0.04 



100 MeV 



cm -2 s _1 MeV" 



(1) 



A radio luminosity similar to the typical values (L r = 10 [erg s 1) of the BL Lac s in the Large 



Area Telescope (LAT) bright AGN sample has been also reported (jAbdo et al.ll2009bl ). but no VHE 
detection above 100 GeV has been reported so far. 

The broad-band SED of NGC 1275 shows similarities with low- frequency-peaked BL Lac (LBL) 



objec ts. VHE 7-ray emission from LBLs is already established, namely from B L Lac (lAlbert et al 



20071 ) and from the two intermediat e-frequency-peaked BL Lac (IBL) 3C 66A (lAcciari et al 



200 



and W Comae ([Acciari et al.ll2008bl ). NGC 1275 is thus a potential VHE source. 



2. VERITAS Observation 



The VERITAS detector is an array of four 12-m diam eter imaging atmospheric-Cherenkov 
telescopes located in southern Arizona (jWeekes et al.l [2002]) . Designed to detect emission from 
astrophysical objects in the energy range from 100 GeV to greater than 30 TeV, VERITAS has 
an energy resolution of ~15% and an angular resolution (68% containment) of ~0.1° per event at 
1 TeV. A source with a flux of 1% of the Crab Nebula flux is detected in less than 50 hours of 
observations, while a 5% Crab Nebula flux source is detected in 2.5 hours. The field of view of 
the VERITAS telescopes is 3.5°. For more details o n the VERITAS instrument and the imaging 
atmospheric-Cherenkov technique, see iHolder et al.l ([2008). 



The Fermi 7-ray space telescope reported to its multiwavelength partners the detection of 
NGC 1275 and a measurement of a hard 7-ray spectrum in the HE range in Fall 2008. VERITAS 
observed the core region of NGC 1275 for about 11 hours between January 15 th and February 26 th 
2009. Data taken with poor weather conditions or technical problems are excluded, resulting in 
7.8 hours of live time. T he average zenith a ngle of the VERITAS observation is 30°. All data were 
taken in "wobble" mode ([Fomin et al.lll994l ) where the telescopes are pointed away from the source 
by 0.5° North/South/East /West to allow for simultaneous background estimation using events from 
the same field of view. 



3. Analysis Methods 



Prior to event selec tion an d bac kgroun d subtraction , all s hower images are calibrated and 
cleaned as described in ICoganl (120061 ) and iDaniel et al.l (|2007l ). Several noise-reducing event- 
selection cuts are made at this point, including rejecting those events where only the two closest- 
spaced telescopes participated in the trigger. Followin g the calibra tion and cleaning of the data, 
the events are parametrized using a moment analysis (jHillaslll985l ). From this moment analysis, 
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scaled parameters are calculat ed and used for the selection of the 7-ray-like events (jAharonian et al 
19971 ; iKrawczvnski et al.ll2006l ). The event selection cuts are optimized a priori on a weak Crab-like 
source (photon index a = 2.5 and 3% Crab Nebula flux level). These selection criteria are termed 
the "standard cuts" . As the VHE spectrum of NGC 1275 is unknown, the possibility that it could 
be considerably softer (photon index a > 4) due to internal (slower bulk motion) or external (dust 
absorption) factors has been considered. A set of modified "soft cuts" are also applied by increasing 
the 9 2 cut (the square of the angular distance of the nominal position of NGC 1275 from the recon- 
structed air shower direction) and decreasing the size cut (the minimum number of photo-electrons 
required in an image). Table [T] shows the set of cuts used in this work. Unless stated otherwise, 
the "standard cuts" are used to generate the results presented in this paper. 



The reflected-region model (jBerge et al.l 120071 ) is used for background subtraction. The total 
number of events in the on-source region is then compared to the total number of events in the 
more numerous off-source regions, scaled by the ratio of the solid angles, to produce a final excess. 



4. Fermi Observation 

The published result from Fermi observation does not overlap with the VERITAS observation 
campaign. Additional data taken during the period from January 15 th to February 26 th 2009 have 
been therefore used in this work. A zenith angle cut at 105° has been applied to eliminate photons 
from the Earth limb. 



The analysis follows the same procedure described in detail in lAbdo et al.1 (|2009al ) . A lower 
energy cut of 200 MeV has been applied to avoid strong systematic uncertainties of the lowest energy 
events. Science Tools version v9rl5p2 and IRFs (Instrumental Response Functions) P6_V3_Diffuse 
(a model of the spatial distribution of photon events calibrated pre-launch) were used throughout 
this paper. The Fermi-LAT data reveal a signal in the energy range between 200 MeV and 6.4 GeV 
of 14.0 a significance at the source location. 



5. Result and discussion 

No significant excess is detected by VERITAS in the NGC 1275 data sample. FigureQ]shows the 
■) 2 plot resulting from the standard cuts analysis. A 99% confide n ce lev el upper limit is calculated 



for both sets of cuts using the method described by iRolke et al.1 (|2005l ) and assuming a Gaussian 
background. The results are reported in Table [T] where the differential flux upper limit at the 
decorrelation energy is shown. We define decorrelation energy the energy that minimizes the 
dependency of the differential upper limit on the assumed photon index. Three different photon 
indices are considered and three "upper limit functions" are hence calculated. The central value of 
the interval where they intersect is taken as decorrelation energy as shown in figure [2j An integral 
flux upper limit above the analysis energy threshold (the energy corresponding to the maximum 
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of the product function of observed spectrum with the collection area, 188 GeV for standard cuts) 
is also shown. The soft-cuts analysis has the effect to lower the analysis energy threshold but also 
to increase the amount of residual background events. Since the decorrelation energy is in both 
analysises above the standard-cuts analysis energy threshold (see table [I]), we discuss hereafter only 
the result from the standard-cuts analysis as it provides a stronger background reduction. 

The energy-binned analysis of the Fermi data taken over the VERITAS observation campaign 
shows a spectrum that can be fitted with the following power law from 200 MeV to 6.4 GeV 



dN 
dE 



(1.74 ±0.45) x 10 



-9 



100 MeV J 



-2.15±0.11 



cm 



- 2 s- 1 MeV" 1 



(2) 



There are some differences from the spectrum presented in lAbdo et al.1 (j2009al ) and summarized 
in formula JT]). The photon index is compatible within the error but the integral flux between 
200 MeV and 6.4 GeV is a factor of 1.37 lower, F^ t = (6.7 ± 0.88) x 10" 8 cm" 2 s~\ compatible 
with the 2008 integral flux within 2.8 standard deviations (F?* t = (9.19 ± 0.50) x 10~ 8 cm" 2 s" 1 ) 



The low er flux comfirms the hint of decreasing trend of the flux that was discussed in lAbdo et al 
(|2009al ). although the difference is statistically weak. In addition, the 2009 spectrum does not 
extend beyond ~6 GeV, because of the overall lower flux and the limited observation time that 
both reduce the statistical significance of the events at higher energy. 

Figure [3] shows the VERITAS differential flux upper limit, calculated at the decorrelation 
energy of the standard-cuts analysis, together with the Fermi spectra. The lower statistical sig- 
nificance of the overall signal i n the Fermi data s imultaneous with the VERITAS observation, 
compared to the measurement in lAbdo et al.l (|2009al ) , is visible in the larger error bars of the spec- 
tral points. The VERITAS upper limit is a factor of five below the extrapolation of the power-law 
fit between 200 MeV and 6.4 GeV of the Fermi data taken during the VERITAS observation cam- 
paign. Considering the proximity of the NGC 1275 galaxy {z = 0.018), negligible absorption by the 
interaction of VHE 7-rays on the extra-galactic background light (EBL) is expected. Given three 
degrees of freedom for the data sample, the VERITAS upper limit implies a power-law fit with a 
normalized x 2 /3 value of 2.97, corresponding to a probability of P{x 2 ) = 0.03. The large value 
of the x 2 probability is given by the large error bars of the spectral p oints due to the rel atively 
short observation time and to the lower flux than the measurement in lAbdo et al.l (|2009a|) , that 



both reduce the overall statistics of events especially at higher energy. Since the photon index 
does not change significantly and the extrapolated flux is still la rger than the VER ITAS flux upper 
limit, the same statistical test has been done on the spectrum in lAbdo et al.l (|2009al ) after applying 
a scaling factor of 1.35 to the spectrum itself. The result gives a normalized x 2 /5 value of 4.57 
with a probability P(x 2 ) = 3.6 x 10~ 4 . We conclude that the upper limit is incompatible with an 
extrapolation of the power law from the Fermi data. 



A deviation from the pure power-law spectrum is suggested. The Fermi data analyzed in this 
work are fitted using a power law with the addition of an exponential cutoff, imposing the constraint 
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given by the upper limit. The VHE result is compatible with this function if the cut-off energy 
is E cut ~ 18 GeV. A normalized x 2 /3 = 1-44 is obtained for the fit. A power law with such an 
exponential cutoff in the IC component may arise from an electron energy distribution that follows 
a power law with a rather sharp cutoff. 

Alternatively, if the energy distribution of the accelerated electrons in the jet follows a power 
law with an exponential cutoff in the form: 

AT (7) oc 7 - p e" 7/7c (3) 

where p is the electron spectral index, 7 is the particle Lorentz factor and j c is the Lorentz factor 
at the cut-off energy, the emitted synchrotron radiation flux and IC component, as a function of 
energy 4>(E), would have a sub-exponential cutoff in the form: 

<p(E) oc ^-(p- 1 )/ 2 e ~V E / E ^ (4) 

This sub-exponential cutoff is a reasonable approximation if the IC scattering proceeds in the 
Thomson regime. Fitting the Fermi-LAT data and including the VERITAS upper limit gives a 
cut-off energy E cu t ~ 120 GeV and x 2 /3 = 1.56. Probabilities of the x 2 quantile with five degrees 
of freedom are 0.23 and 0.20 for the exponential and sub-exponential cut-off result, respectively. 
The fit with a sub-exponential cutoff gives a smaller normalized x 2 value than the fit with the sub- 
exponential cutoff, but the difference is not significant. Both assumptions are therefore compatible 
with Fermi-LAT and VERITAS measurements. 

Another possibility is that the measured power-law steepens to a softer spectral index. This 
behaviour arises if the accelerated electrons below a certain energy escape from their environment 
faster than the typical synchrotron radiation time. Assuming an electron injection in the jet 
following a power law Qij) oc j~ p and a constant and energy- independent leakage of particles per 
unit time, the stationary electron energy distribution may be approximated with the form: 

AT( 7 ) oc 7 - p (l - e^ 7 ) (5) 

where 7 b is the energy at which the energy loss due to the leakage is equal to the synchrotron 
radiating energy loss. The IC component of the SED would assume the following form: 

P{v) = vF(v) oc v^-^l 2 (1 - e-^^^j (6) 

where is the break frequency in the IC component. For i/ < equation ([6]) asymptotically 
follows the power law z/ 1- f)/ 2 , and for v ^> the asymptotic limit becomes v~ p / 2 . Equation ([6]) 
thus represents a smooth transition from a power-law regime below the break frequency to a 
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steeper power-law regime above the break frequency. In this scenario the upper limit obtained by 
the VERITAS measurement implies a break energy E\, ~ 16 GeV. The normalized \ 2 obtained 
from the fit with function ([6]) is x 2 /3 = 1.41. The probability of the x 2 quantile with 3 degrees of 
freedom is 0.24, similar to the probabilities obtained in the exponential- and sub-exponential-cutoff 
scenarios here considered. 

It must be noticed that the result on the cutoff or break energy does not change significantly 
in the soft-cuts analysis since the corresponding upper limit is slightly less constraining due to the 
weaker background reduction. 

6. Conclusions 

VERITAS observed the radio galaxy NGC 1275 (Perseus A) for a total of 7.8 hours good- 
quality live time between January and February 2009. The observation did not result in a detection, 
despite the strong flux and hard spectrum recently reported by Fermi up to 25 GeV during the 2008 
campaign and despite the similar flux during the VERITAS observation campaign. Upper limits at 
99% confidence level are calculated using two analyses optimized for Crab-like and softer spectra, 
respectively. In both cases, the flux upper limit is incompatible with an extrapolation of the power 
law measured by Fermi. Excluding absorption by the EBL, due to the proximity of the galaxy, 
the non-detection by VERITAS might be due to an intrinsic variability of the core. However, the 
Fermi-LAT data during the VERITAS observation campaign show a dim in flux by a factor 1.37 at 
the level of 2.8 a and no significant change in photon index, compared to the 2008 data. A simple 
variability of the flux is, therefore, not sufficient to explain the non detection by VERITAS. The 
possibility of intrinsic spectral curvature at very high energy thus offers a viable explanation of 
the data, but such a spectral variation must generate a strong drop of the flux at an energy not 
much above the highest energy measured by Fermi-L AT . The upper limit in this work suggests an 
exponential cutoff in the Fermi spectrum at an energy E ~ 18 GeV, or E ~ 120 GeV in the case 
of a sub-exponential cutoff. In either case the cutoff is likely above the maximum energy measured 
by Fermi and both assumptions are compatible with the Fermi and VERITAS measurements. The 
possibility that the source spectrum is a broken power law is also considered and is compatible with 
the data with a break energy of about 16 GeV. 

This research was supported by grants from the U.S. Department of Energy, the U.S. National 
Science Foundation and the Smithsonian Institution, by NSERC in Canada, by Science Foundation 
Ireland and by STFC in the UK. 
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Table 1: Selection cuts applied to the VERITAS data for NGC 1275 and the analysis results. For 
each of the two sets of cuts, the nine columns show: the 9 2 cut; the size cut in digital counts 
(1 DC ~ 4.5 photoelectrons); the analysis energy threshold; the decorrelation energy; the tota l 
number of excess events; the corresponding significance according to Formula 17 in lLi &: Mai ([1983); 
the differential upper limit at the decorrelation energy; the same upper limit in Crab Nebula units; 
the integral upper limit above the analysis energy threshold. A total of 9 and 5 background regions 
have been used for the standard and soft cuts analysis respectively. A power-law photon index of 
—2.5 has been assumed in the flux upper limit calculation. It is important to note that the upper 
limits for the standard and soft cuts are highly correlated. 
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Fig. 1.— 9 2 plot for VERITAS observation of NGC 1275. The points with error bars represent the 
ON-source data sample and the filled area is the background. 
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Fig. 2. — Decorrelation energy for the standard-cuts analysis on VERITAS data. Three different 
photon indices, —2.17, —3.0 and —4.0, are assumed for the power law used to calculate the flux 
upper limit. The central value of the interval where the three normalized power laws intersect is 
taken as decorrelation energy. The dependency of the differential flux upper limit on the assumed 
photon index is strongly reduced at the decorrelation energy. Assuming a photon index —2.5, the 
effect of the photon index change within the values here considered at the decorrelation energy is 
about 15%. 
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• Fermi-LAT data (this work) 
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Fig. 3. — NGC 1275 spectrum and the VERITAS upper limit on the differential flux at the decor- 
relation energy 338 GeV (standard cuts). The solid circles with error bars are the measurement by 
the Fermi 7-ray space telescope during the VERITAS o b servati on campaign. Empty circles with 
error bars are the measurement presented in lAbdo et al.l (|2009al ) from the energy-binned analysis. 
The solid line is the power-law fit to the Fermi data. The dashed line is the extrapolation of the 
power-law. The dotted-dashed line is the fit of a power law with an exponential cutoff at 18 GeV. 
The double-dotted dashed line is the fit of a power law with a sub-exponential cutoff at 120 GeV 
and the dotted line is the smooth broken power law fit of a break energy at 16 GeV. All fits are 
done on the Fermi data analyzed in this work (solid circles). 



